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1 Introduction 



The totally inclusive cross section of deep-inelastic lepton-proton scattering (DIS) depends 
on the square s of the centre-of-mass energy, Bjorken's variable x = Q'^/{2pq), and the 
inelasticity variable y = Q'^/{xs), where p and q are the four-momenta of the proton and 
the virtual photon, respectively, and = —q^ > 0. The doubly differential cross section 
is parameterized in terms of the structure function F2 and the longitudinal structure 
function Fl, as 

^""""^ll + (1 - ymix, Q') - y'F.ix, Q')}, (1) 



dx dy xQ^ 

where a is Sommerfeld's fine-structure constant. At small values of a;, F^ becomes non- 
negligible and its contribution should be properly taken into account when the F2 is 
extracted from the measured cross section. The same is true also for the contributions F2 
and Fl of F2 and Fl due to the heavy quarks i = c,b. 

Recently, the HI P El |3| and ZEUS ^ El E] Collaborations at HERA presented new 
data on F2 and F2 . At small x values, of order 10~^, F2 was found to be around 25% of F2, 
which is considerably larger than what was observed by the European Muon Collaboration 
(EMC) at CERN [7] at larger x values, where it was only around 1% of F2. Extensive 
theoretical analyses in recent years have generally served to establish that the F^ data can 
be described through the perturbative generation of charm within QCD (see, for example, 
the review in Ref. [S] and references cited therein). 

In the framework of Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) dynamics 
[9], there are two basic methods to study heavy- flavour physics. One of them [10] is based 
on the massless evolution of parton distributions and the other one on the photon-gluon 
fusion (PGF) process [12j. There are also some interpolating schemes (see Ref. [13] and 
references cited therein). The present HERA data on ^2^^ [H El El IH ISl E] are in good 
agreement with the modern theoretical predictions. 

In earlier HERA analyses [H H] , i^f and F^ were taken to be zero for simplicity. Four 
years ago, the situation changed: in the ZEUS paper the F£ contribution at next-to- 
leading order (NLO) was subtracted from the data; in Refs. [21 E], the HI Collaboration 
introduced the reduced cross sections 

for i = c,b and thus extracted at NLO by fitting their data. Very recently, a sim- 
ilar analysis, but for the doubly differential cross section rf^a'Y {dx dy) itself, has been 
performed by the ZEUS Collaboration 

In this letter, we present a compact formula for the ratio Ri = FI/F2, which greatly 
simplifies the extraction of F2 from measurements of d'^a^^/{dxdy). 
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2 Master formula 



We now derive our master formula for Ri{x, Q"^) appropriate for small values of x, which 
has the advantage of being independent of the parton distribution functions (PDFs) 
fa{x,Q^), with parton label a = g,q,q, where q generically denotes the light-quark 
flavours. In the small-x range, where only the gluon and quark-singlet contributions 
matter, while the non-singlet contributions are negligibly small, we hav^ 

Fl{x,Q')= 5^ Ci>,Q2)®a;/^(x,Q^), (3) 

a=g,q,ql=+,- 

where / = ± labels the usual + and — linear combinations of the gluon and quark- 
singlet contributions, Cl ^{x, Q"^) are the DIS coefficient functions, which can be calculated 
perturbatively in the parton model of QCD, /i is the renormalization scale appearing in 
the strong-coupling constant ^^(/i), and the symbol ® denotes convolution according to 
the usual prescription, /(x) ® g{x) = J^{dy/y)f{y)g{x/y). Massive kinematics requires 
that C{^^ = for X > 6j = l/(l+4aj), where a, = mf/Q'^. We take rrii to be the solution of 
^iiiT^i) = fni, where mi{n) is defined in the modified minimal-subtraction (MS) scheme. 
Exploiting the small-x asymptotic behaviour of fj^i^XjQ"^) [12], 

(4) 

Eq. ([3]) can be rewritten as 

FUx,Q')^ J2 5^ M^1 + 5,,Q>/^(x,Q2), (5) 

a.=9,q,ql=+,- 

where ^ 

Ml,{n,Q')= r dxx--'Cl,{x,Q^) (6) 

^0 

is the Mellin transform, which is to be analytically continued from integer values n to real 
values 1 + 5i. 

As demonstrated in Ref. [13], HERA data support the modified Bessel-like behav- 
ior of PDFs at low X values predicted in the framework of the so-called generalized 
double-asymptotic scaling regime. In this approach, one has M^^(1,Q^) = M^^(1,Q^) 
if M^^(n, Q^) are devoid of singularities in the limit 5/ — > 0, as in our case. Defining 
Mfc,,(l,Q2) = M±„(l,g2) and using /„(x,Q2) = Y.i=± fii^.Q^), Eq. © may be simpli- 
fied to become 

Fl{x,Q^)^ J2 M,,ail,Q')xfaix,Q'). (7) 
a-=g,q,q 

A further simplification is obtained by neglecting the contributions due to incoming light 
quarks and antiquarks in Eq. ([7]), which is justified because they vanish at LO and are 

^Here and in the following, we suppress the variables fi and rrii in the argument lists of the structure 
and coefficient functions for the ease of notation. 
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numerically suppressed at NLO for small values of x. One is thus left with the contribution 
due to PGF [E], 

F^(x,Q2) ~ Mu,g{l,Q^)xfg{x,Q^). (8) 
In fact, the non-perturbative input fg{x,Q'^) does cancels in the ratio 



Ri{x,Q 



(9) 



which is very useful for practical applications. Through NLO, Mfcg(l,Q^) exhibits the 
structure 



M,,,(l, Q^) = efai^) <» M^^Vl, a.) + a(/i) 



Ml]l{l,a.)+Mi'>(l,a,) 



1 

X In — ^ 

mf 



(10) 



where Cj is the fractional electric charge of heavy quark i and a(/i) = as(/i)/(47r) is the 
couplant. Inserting Eq. ffTOj) into Eq. ([9]), we arrive at our master formula 



MS(l,a.) 


+ a(/i) 






ln(/i^/mf) 


Mi°)(l,a,) 


+ a(/i) 




+ Mg(l,a.) 


ln(yU^/mf ) 



+ 0(a2). 



We observe that the right-hand side of Eq. fill I) is independent of x, a remarkable feature 
that is automatically exposed by our procedure. In the next two sections, we present 
compact analytic results for the LO (j = 0) and NLO (j = 1,2) coefficients M^"'^(l, a^), 
respectively. 



3 LO results 

The LO coefficient functions of PGF can be obtained from the QED case [15j by adjusting 
coupling constants and colour factors, and they read |16l [17] 

Ci^l{x, a) = -2x{[l - 4x(2 - a)(l - x)]p - [1 - 2a;(l - 2a) 

+ 2x^(1- 6a -Aa^)]L{P)}, 
d^lix, a) = 8x^[{l - x)(3 - 2axL{(3)], (12) 

where 



m^u.l±l. (13) 
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Using the auxiliary formulas 



^0 



1 — 2a J (a), if m = 

a;-/5= <( |[1 -2a-4a(l + 3a)J(a)], if m = 1 , (14) 



I 
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[(1 + 3a)(l + 10a) - 6a(l + 6a + lOa^) J(a)], if m = 2 



„o I J (a), if m = 

/ a;'"L(/3) = i -|[1- (1 + 2a) J(a)], if m = 1 , (15) 

' -f [3(l + 2a)-2(l + 4a + 6a2)J(a)],if m = 2 



where 



J{a) = -Vb\nt, t = - — (16) 

1 + v6 



we perform the Mellin transformation in Eq. ([6]) to find 

m£)(1, a) = ^[1 + 2(1- a) J(a)], 



m{,°J(1, a) = h[l + 6a - 4a(l + 3a) J(a)]. (17) 
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At LO, the small-x approximation formula thus reads 

1 + 6aj — 4aj(l + 3aj) J(a, 



1 + 2(1 -ai)J(a,) ■ ^ ' 

4 NLO results 

The NLO coefficient functions of PGF are rather lengthy and not published in print; they 
are only available as computer codes [18] . For the purpose of this letter, it is sufficient to 
work in the high-energy regime, defined by a^ -C 1, where they assume the compact form 

m 

Cg(x,a)=/?/?£(l,a), (19) 



with 



i?(^)(l, a) = -Ca[5 + (13 - 10a) J(a) + 6(1 - a)/(a)], 
9 

R'-ini, a) = -^CAb{l - 12a - [3 + 4a(l - 6a)] J(a) + 12a(l + 3a)/(a)}, 



'>9^ ' ' g 

4'J(l,a) = -4C^M(°)(l,a), (20) 
where Ca = N for the colour gauge group SU(N), J(a) is defined by Eq. ( |T6l) . and 

1 



Ua) = -Vb 



C(2) + I In^ t - ln(a6) Int + 2 Li2(-t) 



(21) 



Here, ({2) = 7r^/6 and Li2(x) = — J^{dy/y) \n{l — xy) is the dilogarithmic function. Using 
Eq. (IT^ for m = 0, we find the Mellin transform of Eq. (fT9l) to be 

Mi^)(l,a) = [l-2aJ(a)]i?l;'J(l,a). (22) 
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Table 1: Values of F^lx, Q"^) extracted from the HI measurements of a'^'^ at low [3] and 
high P] values of (in GeV^) at various values of x (in units of 10~^) using our approach 
at NLO for /i^ = ^Q"^ with ^ = 1, 100. The LO results agree with the NLO results for 
^ = 1 within the accuracy of this table. For comparison, also the results determined in 
Refs. [21 [3] are quoted. 





X 


HI 




= 100 g2 


12 


0.197 


0.435 ± 0.078 


0.433 


0.432 


12 


0.800 


0.186 ±0.024 


0.185 


0.185 


25 


0.500 


0.331 ±0.043 


0.329 


0.329 


25 


2.000 


0.212 ±0.021 


0.212 


0.212 


60 


2.000 


0.369 ± 0.040 


0.368 


0.368 


60 


5.000 


0.201 ±0.024 


0.200 


0.200 


200 


0.500 


0.202 ± 0.046 


0.201 


0.201 


200 


1.300 


0.131 ±0.032 


0.130 


0.130 


650 


1.300 


0.213 ±0.057 


0.212 


0.213 


650 


3.200 


0.092 ±0.028 


0.091 


0.091 



5 Results 

We are now in a position to explore the phenomenological implications of our results. 
As for our input parameters, we choose rric = 1.25 GeV and mj, = 4.2 GeV. While the 
LO result for Ri in Eq. (fT8l) is independent of the unphysical mass scale /i, the NLO 
formula (|TT1) does depend on it, due to an incomplete compensation of the /i dependence 
of a(/i) by the terms proportional to ln(/x^/(5^), the residual dependence being formally 
beyond NLO. In order to estimate the theoretical uncertainty resulting from this, we put 
yU^ = C,Q^ and vary ^. Besides our default choice ^ = 1, we also consider the extreme 
choice ^ = 100, which is motivated by the observation that NLO corrections are usually 
large and negative at small x values [20]. A large ^ value is also advocated in Ref. [21] . 
where the choice ^ = 1/x", with 0.5 < a < 1, is proposed. 

We now extract F2{x,Q'^) {i = c,b) from the HI measurements of the reduced cross 
sections in Eq. at low {12 < < 60 GeV^) [5] and high {Q^ > 150 GeV^) [2] values 
of using the LO and NLO results for Ri derived in Sections [3] and HJ respectively. 
Our NLO results for /x^ = ^Q^ with ^ = 1, 100 are presented for i = c,b in Tables [T] and 
|2l respectively, where they are compared with the values determined by HI. We refrain 
from showing our results for other popular choices, such as /i^ = 4m^, + 4mf because 
they are very similar. We observe that the theoretical uncertainty related to the freedom 
in the choice of /i is negligibly small and find good agreement with the results obtained 
by the HI Collaboration using a more accurate, but rather cumbersome procedure [21 E]. 
The experimental data from the ZEUS Collaboration [6] do not allow for such an analysis 
because they do not come in the form of Eq. ([2]) . 

In order to assess the significance of and the theoretical uncertainty in the NLO cor- 
rections to Ri, we show in Fig. [Tjthe dependences of Rc, Rb, and Rt evaluated at LO 
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Table 2: Values of -F2^(x, Q"^) extracted from the HI measurements of a^^ at low [3] and 
high [2] values of (in GeV^) at various values of x (in units of 10~^) using our approach 
at NLO for /x^ = ^j^j^ ^ = 1, 100. The LO results agree with the NLO results for 
^ = 1 within the accuracy of this table. For comparison, also the results determined in 
Refs. [21 [3] are quoted. 





X 


HI 




/i^ = 100 Q2 


12 


0.197 


0.0045 ±0.0027 


0.0047 


0.0046 


12 


0.800 


0.0048 ±0.0022 


0.0048 


0.0048 


25 


0.500 


0.0123 ±0.0038 


0.0124 


0.0124 


25 


2.000 


0.0061 ±0.0024 


0.0061 


0.0061 


60 


2.000 


0.0190 ±0.0055 


0.0190 


0.0190 


60 


5.000 


0.0130 ±0.0047 


0.0130 


0.0130 


200 


0.500 


0.0413 ±0.0128 


0.0400 


0.0400 


200 


1.300 


0.0214 ±0.0079 


0.0212 


0.0212 


650 


1.300 


0.0243 ±0.0124 


0.0238 


0.0238 


650 


3.200 


0.0125 ±0.0055 


0.0125 


0.0125 



from Eq. ( |T8l) and at NLO from Eq. (fTTi) with /i^ = Amf,Q'^ ± 4mf. We observe from 
Fig. [1] that the NLO predictions are rather stable under scale variations and practically 
coincide with the LO ones in the lower regime. On the other hand, for ^ 4m^, the 
NLO predictions overshoot the LO ones and exhibit an appreciable scale dependence. We 
encounter the notion that the fixed-flavour-number scheme used here for convenience is 
bound to break down in the large-Q^ regime due to unresummed large logarithms of the 
form ln{Q'^/mf). In our case, such logarithms do appear linearly at LO and quadratically 
at NLO. In the standard massless factorization, such terms are responsible for the 
evolution of the PDFs and do not contribute to the coefficient functions. In fact, in the 
variable-flavour-number scheme, they are MS-subtracted from the coefficient functions 
and absorbed into the evolution of the PDFs. Thereafter, the asymptotic large-Q^ 
dependences of Ri at NLO should be proportional to ^^(Q^) and thus decreasing. This is 
familiar from the Callan-Gross ratio R = Fl/{F2 — Fl), as may be seen from its (x, Q^) 
parameterizations in Ref. [22]. Fortunately, this large-Q^ problem does not affect our re- 
sults in Tables [1] and [2] because the bulk of the HI data is located in the range of moderate 
values. Furthermore, Ri enters Eq. ([2]) with the suppression factor y'^/[l ± (1 — vY]- 
The ratio Rc was previously studied in the framework of the /cj-factorization approach 
[T7] and found to weakly depend on the choice of unintegrated gluon PDF and to be 
approximately x independent in the small-x regime (see Fig. 8 in Ref. [T7]). Both features 
are inherent in our approach, as may be seen at one glance from Eq. f fTT]) . The prediction 
for Rc from Ref. [T7], which is included in Fig.[T]for comparison, agrees well with our results 
in the lower range, but it continues to rise with Q^, while our results reach maxima, 
beyond which they fall. In fact, the fc^-factorization approach is likely to overestimate Rc 
for ^ 4m^, due to the unresummed large logarithms of the form hi{Q'^/mj) discussed 
above. 
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Figure 1: Rc, Rb, and Rt evaluated as functions of at LO from Eq. ( fTSi) (dot-dashed 
lines) and at NLO from Eq. ( |TT|) with //^ = imf (dashed lines) and /x^ = + 4mf (solid 
lines). For comparison, the prediction for Rc in the /ct-factorization approach (dot-dot- 
dashed line) [17] is also shown. 
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6 Conclusions 



In this letter, we derived a compact formula for the ratio Ri = FI/F2 of the heavy- 
flavour contributions to the proton structure functions F2 and Fl valid through NLO at 
small values of Bjorken's x variable. We demonstrated the usefulness of this formula by 
extracting F2 and F| from the doubly differential cross section of DIS recently measured by 
the HI Collaboration [21 [3] at HERA. Our results agree with those extracted in Refs. [21 [3] 
well within errors. In the range probed by the HI data, our NLO predictions agree 
very well with the LO ones and are rather stable under scale variations. Since we worked 
in the fixed-fiavour-number scheme, our results are bound to break down for ^ 4m^, 
which manifests itself by appreciable QCD correction factors and scale dependences. As is 
well known, this problem is conveniently solved by adopting the variable-fiavour-number 
scheme, which we leave for future work. Our approach also simply explains the feeble 
dependence of Ri on x and the details of the PDFs in the small-x regime. 
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